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We review the infrared behavior of interacting bosons at zero temperature. After a brief 
discussion of the Bogoliubov approximation and the breakdown of perturbation theory 
due to infrared divergences, we show how the non-perturbative renormalization group 
enables to obtain the exact infrared behavior of the correlation functions. 
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1. Introduction 

Many of the predictions of the Bogoliubov theory of superfluidi ty-^ have been con- 
firmed experimentally, in particular in ultracold atomic gases.L^ Nevertheless a 
clear understanding of the infrared behavior of interacting bosons at zero temper- 
ature has remained a challenging theoretical issue for a long time. Early attempts 
to go beyond the Bogoliubov theory have revealed a singular perturbation theory 
plagued by infrared diver gence s due to the presence of the Bose-Einstein condensate 
and the Goldstone modeSEEIln the 1970s, Nepomnyashchii and Nepomnyashchii 
proved that the anomalous self-energy vanishes at zero frequency and momentum 
in dimension d < 3. This exact result shows that the Bogoliubov approximation, 
where the linear spectrum and the superfluidity rely on a finite val ue o f the anoma- 
lous self-energy, breaks down at low energy. As realized latter onJ^iSI the singular 
perturbation theory is a direct consequence of the coupling between transverse and 
longitudinal fluctuations and reflects the divergence of the longitudinal susceptibil- 
ity - a general phenomenon in systems with a continuous broken symmetryP^ 

In this paper, we review the infrared behavior of interacting bosons. A more 
detailed discussion together with a comparison to the classical O(N) model can 
be found in Ref. In Sec. [3J we briefly review the Bogoliubov theory and 
the appearance of infrared divergences in perturbation theory. We introduce the 
Ginzburg momentum scale pa signaling the breakdown of the Bogoliubov ap- 
proximation. In Sec. [31 we discuss the non-perturbative renormalization group 
(NPRG) and show how it enables to obtain the exact infrared behavior of the 
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normal and anomalous single-particl e propagators without encountering infrared 
divcrgcnccc .^ 13114115 ' 1611 ^ 118 ' 191201 



2. Perturbation theory and breakdown of the Bogoliubov 
approximation 

We consider interacting bosons at zero temperature with the (Euclidean) action 
S = J dx ^*(d T -a-^)^ + |(^» 2 , (1) 

where ip( x ) is a bosonic (complex) field, x — (r, r), and J dx = f£ dr J d d r. r <E 
[0, ft] is an imaginary time, ft — > oo the inverse temperature, and fi denotes the 
chemical potential. The interaction is assumed to be local in space and the model 
is regularized by a momentum cutoff A. We consider a space dimension d > 1. 
Introducing the two-component field 

*(p) = (rl-p)) ' ¥{p) = (^(p).^(-p)) (2) 

(with p — (p, iu) and uj a Matsubara frequency), the one-particle (connected) 
propagator becomes a 2 x 2 matrix whose inverse in Fourier space is given by 



(3) 



J + (J, - e p - E n (p) -San(p) 

- S anb) -ZW + /X - £ p - E a (-p) 

where S n and S an are the normal and anomalous self-energies, respectively, and 
e p = p 2 /2m. If we choose the order parameter (ip(x)) — ^Jno to be real (with no 
the condensate density), then the anomalous self-energy E an (p) is real. 



2.1. Bogoliubov approximation 

The Bogoliubov approximation is a Gaussian fluctuation theory about the saddle 
point solution tp{x) = i/n~o = \J u/g. It is equivalent to a zero-loop calculation of 
the self-energies! 21 ! 22 ! 

^ n °Hp) = 2gn , ^(p) = gn . (4) 
This yields the (connected) propagators 

G$»<p) = -WpW(p)) c - ~ iu -^~ 9no , 

CO + £/ p 

G^{p) = -<V(P)^(-P))c = ( 5 ) 

where i? p = [e p (e p + 2gno)] 1 / 2 is the Bogoliubov quasi-particle excitation energy. 
When |p| is larger than the healing momentum p c = (2gmno) 1 ' 2 , the spectrum 
E p ~ e p + griQ is particle-like, whereas it becomes sound-like for |p| -C p c with a 
velocity c = \Jg~noJm. In the weak-coupling limit, no — n (n is the mean boson 
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density) and p c can equivalently be defined as p c — (2gmn) 1 ^ 2 . In the Bogoliubov 
approximation, the occurrence of a linear spectrum at low energy (which implies 
superfluidity according to Landau's criterion) is due to E an (p = 0) being nonzero. 



2.2. Infrared divergences and the Ginzburg scale 

The lowest-order (one-loop) correction E^ to the Bogoliubov result S'- -' is diver- 
gent for d < 3. Retaining only the divergent part, we obtain 

-i m 3 / ,,2\(<*-3)/2 



if d < 3 and 



4 3 / 2 \ \ a ~ °) l * 

E«(p) * E«(p) * ~^A d+1 ( P 2 + ^ j (6) 



E«(p) * E«(p) * -^Uln ( ^ ) (7. 



if <i = 3, where 



^M2)ifd<4, 



A d = < sin(7rd/2) r(d-l) 11 " ^ ^ ,gx 

if d = 4. 

We can estimate the characteristic (Ginzburg) momentum scale pa below which 
the Bogoliubov approximation breaks down from the condition lE^^p)! ~ T,^\p) 
or |E a n (p)| — Ein (p) for |p| = p G and |w| = cp G , 

f {A d+1 gmp c ) l l^ if d< 3, 

This result can be rewritten as 

U c (A d+1 5' i / 2 ) 1 /(3-rf)if d <3, 
jp c exp(-^i^) if d = 3, 

where 

~g = gmn^ ~ (J^) 2 (11) 

is the dimensionless coupling constant obtained by comparing the mean interaction 
energy per particle gfi to the ty pica l kinetic energy 1 jmf 2 where f ~ is the 

mean distance between particles!^ 1 A superfluid is weakly correlated if ^ <gc 1, i.e. 
Pg *C p c <C fi 1 /'' (the characteristic momentum scale n 1 /" 1 does however not play 
any role in the weak-coupling limit). In this case, the Bogoliubov theory applies to 
a large part of the spectrum where the dispersion is linear (i.e. |p| < p c ) and breaks 
down only at very small momenta |p| < pc *C p c - When the dimensionless coupling 
g becomes of order unity, the three characteristic momentum scales pc ~ Pc ~ n. 1 ^ 
become of the same order. The momentum range [pcPc] where the linear spectrum 
can be described by the Bogoliubov theory is then suppressed. We expect the strong- 
coupling reg ime g ^> 1 to be governed by a single characteristic momentum scale, 
namely n 1/,£Z P31 



4 N. Dupuis 



2.3. Vanishing of the anomalous self-energy 

The exact values of £ n (p = 0) and £ an (p = 0) can be obtained using the U(l) 
symmetry of the action, i.e. the invariance under the field transformation if)(x) — > 
e ie ip(x) and ip*( x ) ~> e~ l6 _ib* (x) . On the one hand, the self-energies satisfy the 
Hugenholtz-Pines theorem,^ 

£ n (p = 0) - S an (p = 0) = a. (12) 

On the other hand, the anomalous self-energy vanishes, 

£ an Cp = 0)=0. (13) 

The last result was first proven by Nepomnyashchii and Nepomnyashchii JZEQUI] 
shows that the Bogoliubov theory, where the linear spectrum and the superfluidity 
rely on a finite value of the anomalous self-energy, breaks down at low energy in 
agreement with the conclusions drawn from perturbation theory (Sec. 12 .2[) . 



3. The non-perturbative RG 

The NPRG enables to circumvent the difficulties o f perturbation theory and de- 
rive the correlation functions in the low-energy li m itP^^ 3 | 16 | 17 | 18 | 14 | 15 | 19 | 2 0] The 
strategy of the NPRG is to build a family of theories indexed by a momentum scale 
k such that fluctuations are s moothly taken into account as k is lowered from the 
microscopic scale A down to oPSESl This is achieved by adding to the action (jTJ) 
an infrared regulator term 

As k [r,i>] = J2^(p) R k(p)^(p)- (14) 
p 

The main quantity of interest is the so-called average effective action 

r fc [0* ,<!>] = - in z k [j* ,J] + J2 1 J * (p)^(p) + c - c -l ~ AS * W > 0. ( 15 ) 

p 

defined as a modified Legendre transform of — In Zk [J* , J] which includes the sub- 
traction of ASk [(/)*, (f>]. J denotes a complex external source that couples linearly to 
the boson field ip and 4>[x) — (ip(x)) is the superfluid order parameter. The cutoff 
function Rk is chosen such that at the microscopic scale A it suppresses all fluctu- 
ations, so that the mean-field approximation Ta[4>*,<P] — S [(/)*, </)} becomes exact. 
The effective action of the original model ([T]) is given by Tk=o provided that Rk=o 
vanishes. For a generic value of fc, the cutoff function Rk(p) suppresses fluctuations 
with momentum |p| < fc and frequency \oj\ < ck but leaves those with |p|, \oj\/c > k 
unaffected (c = ct is the velocity of the Goldstone mode) . The dependence of the 
average effective action on k is given by Wetterich's equatiorPSl 

5 t r fc [0*,0] = iTr|7? fc (r^[r,0] +J R fc ) _1 |, (16) 
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(2) 

where t — ln(fc/A) and Rk = d t Rk- 1^ [4>*,(j)] denotes the second-order functional 
derivative of [</>]. In Fourier space, the trace involves a sum over momenta and 
frequencies as well as the internal index of the <fi field. 



3.1. Derivative expansion and infrared behavior 

Because of the regulator term AS kl the vertices T k (p%, ■ ■ ■ ,p n ) are smooth func- 
tions of momenta and frequencies and can be expanded in powers of pf/k 2 and 
Lo 2 /c 2 k 2 . Thus if we are interested only in the low- energy properties, we can use a 
derivative expansion of the average effective actionf2S126| j n f n ow i n g we con- 
sider the ansat 

\z c , k d T -V A . k d 2 T ZA > k - 



dx 



2m 



2 V 



(17) 



where n = \4>\ 2 - ^o,fe denotes the condensate density in the equilibrium state. We 
have introduced a second-order time derivative term. Although not present in the 
initial aver age effe ctive action T\, we shall see that this term plays a crucial role 
when d < 

In a broken U(l) symmetry state with real order parameter <p — ^/uq, the normal 
and anomalous self-energies are given by 

Sfc, n (p) = M + Vam^ 2 + (1 - Z c ,k)iu - (1 - Z A ,k)e P + Afcn 0j fc, (18) 

Sfe,an(p) = A/-no,fe. 



These expressions imply the existence of a sound mode with velocity 

/ Z A . k /2m \ 

Ck 



V A ,k + Zq k /2X k n 0y k 



(19) 



At the initial stage of the flow, Za.a — Zc,a — 1, V a ,a — 0, Aa = g and no,A = m/.9, 
which reproduces the results of the Bogoliubov approximation. A crucial property 
of the RG flow is that 



A* 



,3-o 



(20) 



vanishes with k when d < 3. Eq. (120|) follows from the numerical solution of the 
RG equations, but can also be anticipated from the expected singular behavior of 
the longitudinal propagator!!!! 

The parameters ^ M.fc, Zc.k and Va,u can be related to thermodynamic quantities 
using Ward identities p | 13 | 18 | 



n s ,k = ZA,kno yk = n k , 
1 d 2 U k 



V A 



Zc,k 



2n ,fc 9/i 2 
d 2 U k 



(21) 



dn. 



Q.k 



6 N. Dupuis 




Fig. 1. (Color online) (Left panel) A*,, Zq j. and Va k vs m (PG/k) where pq = y (gm) 3 n/Air for 
n = 0.01, 2mg = 0.1 and d = 2 [ln(pQ/p c ) ~ —5.87]. The inset shows vs 2mg obtained from 
the criterion Va, pg = VA,k=o/2 [the Green solid line is a fit to pq ~ (2mg) 3 ' 2 ]. (Right panel) 
Condensate density no fe, superfluid density n s /, and Goldstone mode velocity c^- vs \n(pQ/k). 



where fik is the mean boson density and n s ^ the superfluid density. Here we con- 
sider the effective potential Uk as a function of the two independent variables n 
and \x. The first of equations (|21[) states that in a Galilean invariant superfluid at 
zero temperature, the superfluid density is given by the full density of the fluid.^ 
Equations (|2"Tj) also imply that the Goldstone mode velocity Ck coincides with the 
macroscopic sound 

velocity jsn3HHi ie . 

Mk = (22) 

Since thermodynamic quantities, including the condensate "compressibility" 
dno^k/d/j, should remain finite in the limit k — ¥ 0, we deduce from pip that 
Zc,k ~ ^fc ~ k 3 ~ d vanishes in the infrared limit, and 



lim c k = lim . (23) 

Both Zj± k = nk/no,k and the macroscopic sound velocity Ck being finite at k = 0, 
VA.fc (which vanishes in the Bogoliubov approximation) takes a non-zero value when 
k -)• 0. 

The suppression of i?c,fcj together with a finite value of Va,/c=o shows that the 
average effective action (|17|) exhibits a "relativistic" invariance in the infrared limit 
and therefore becomes equivalent to that of the classical 0(2) model in dimensions 
d + 1. In the ordered phase, the coupling constant of this model vanishes as ~ 
fc4-(<m)jn] w hich agrees with (f2"U|). For k — >• 0, the existence of a linear spectrum is 
due to the relativistic form of the average effective action (rather than a non-zero 
value of \kn>o,k as in the Bogoliubov approximation). 

To obtain the fc = limit of the propagators (at fixed p), one should in principle 



stop the flow when k ~ y/p 2 + w 2 /c 2 l^ Since thermodynamic quantities are not 
expected to flow in the infrared limit, they can be approximated by their k = 
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values. Making use of the Ward identities (|21j) . we deduce the exact infrared 
behavior of the normal and anomalous propagators (at k — 0)p^ 

n mc 2 1 mc 2 dn iu 1 



Ca(p) = nomc2 — L - I G|| b), (24) 



where 



G ll (P) - 2n C( W 2 + c 2 p 2)(3-d)/2 ( 25 ) 

is the propagator of the longitudinal fluctuations. The constant C follows from 
the replacement A fe — > C(u; 2 + c 2 p 2 )( 3 ~ d )/ 2 . The leading terms in (j2~4")l agree with 
the results of Gavoret and NozieresPThe contribution of the diverging longitud inal 
correlation function was first identified by Nepomnyashchii and Nepomnyashchii 



3.2. RG flows 

The conclusions of the preceding section can be obtained more rigorously from the 
RG equation (TTB)) satisfied by the average effective action. The flow of Afc, Zc,k and 
Va,u is shown in Fig. [1] for a two-dimensional system in the weak-coupling limit. 
We clearly see that the Bogoliubov approximation breaks down at a characteristic 
momentum scale pa ~ y (gm) 3 n. In the Goldstone regime k <C pg, w e find that 
both Afe and Zc,k vanish linearly with k in agreement with the conclusions of 
Sec. 13.11 Furthermore, Va,u takes a finite value in the limit k — > in agreement 
with the limiting value (|23[) of the Goldstone mode velocity. Figure Q] also shows the 
behavior of the condensate density no,fc, the superfluid density n St k — ZA,kfio,k an d 
the velocity c^. Since ZA,k=o — 1-004, the mean boson density hk = n s- k is nearly 
equal to the condensate density Uo,fc. Apart from a slight variation at the beginning 
of the flow, no,fc, n s _k = ZA,kno,k and Ck do not change with k. In particular, they 
are not sensitive to the Ginzburg scale pq. This result is quite remarkable for the 
Goldstone mode velocity Cfc, whose expression (|19j) involves the parameters Afc, Zc,k 
and Va^i which all strongly vary when k ~ pq. These findings are a nice illustration 
of the fact that the diverg ence of the longitudinal susceptibility does not affect local 
gauge invariant quantities {EEH1 



4. Conclusion 

Interacting bosons at zero temperature are characterized by two momentum scales: 
the healing (or hydrodynamic) scale p c and the Ginzburg scale pa- Pg sets the 
scale at which the Bogoliubov approximation breaks down. For momenta |p| <C p c , 
it is possible to use a hydrodynamic description in terms of density and phase 
variables. This description allows one to compute the correlation functions without 
encountering infrared divergences! ^"^ In this paper, we have reviewed another 
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approach, based on the NPRG, which enables to describe the system at all energy 
scales and yields the exact infrared behavior of the single-particle propagator. A nice 
feature of the NPRG is that it can be used to study models of strongly-correlated 
bosons such as the Bose-Hubbard modelPS 
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